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Abstract 

We introduce a novel type of phase diagram for black holes and black strings on cylinders. 
The phase diagram involves a new asymptotic quantity called the relative binding energy. 
We plot the uniform string and the non-uniform string solutions in this new phase diagram 
using data of Wiseman. Intersection rules for branches of solutions in the phase diagram 
are deduced from a new Smarr formula that we derive. 



1 Introduction 



Neutral and static black holes and black strings on cylinders W^^^ x have proven to 
have a very rich phase structure, a structure which is still largely unknown. Gregory and 
Laflamme |^ H] discovered that uniform black strings on cylinders, i.e. strings that are 
symmetric around the cylinder, are unstable to linear perturbation when the mass of the 
string is below a certain critical mass. This was interpreted to mean that a light uniform 
string decays to a black hole on a cylinder since that has higher entropy. 

However, Horowitz and Maeda |H] argued that such a transition has an intermediate 
step: The light uniform string decays to a non- uniform string, which then eventually 
decays to a black hole. This prompted a search for this missing link, and a branch of 
non-uniform string solutions was found in [1J[3|.^ Furthermore, in |7| the classical decay 
of the Gregory-Laflamme instability was studied explicitly. 

A related question is what happens to a black hole on a cylinder when one increases the 
mass. Obviously the black hole must have a phase transition into a black string when the 
black hole horizon becomes so big that it cannot fit on the cylinder. Several suggestions for 
the phase structure of the Gregory-Laflamme instability, the black hole/string transitions 
and the uniform/non- uniform string transitions have been put forward [HI IH IH} IH} TH ] 
1121 113j but it is unclear which of these, if any, are correct.^ 

In this paper we give new tools to understand the phase structure of neutral and static 
black holes and strings on cylinders. As part of this, we introduce a new type of phase 
diagram for black holes and strings on cylinders. The phase diagram has two parameters, 
the mass M and a new asymptotic quantity called the relative binding energy n. We 
explain why this phase diagram is natural to use from several points of view, and how it 
is useful for study of thermodynamics. We plot moreover the uniform string branch and 
the non-uniform string branch of Wiseman in the phase diagram. 

In Section we consider the linearized Einstein equations for a localized distribu- 
tion of static and neutral matter on a cylinder M'^^^ x S^. In particular we introduce 
binding energy along the periodic direction. From this we define the total mass M and 
the so-called relative binding energy n. We furthermore define the Fourier modes of the 
mass-distribution. We explain that these quantities completely characterize the energy 
momentum tensor. Finally, we show how to measure these quantities from the metric 
independent of the gauge. In particular, this enables us to find M and n for any black 
hole and string solution on the cylinder from the asymptotics of the metric. 

In Section|31we introduce a new phase diagram for static and neutral black holes/strings 
with the mass M and relative binding energy n as the variables. We plot the three known 
branches of solutions in this phase diagram. This includes the black hole branch and the 
uniform string branch, which we explain has n = l/{d — 2). We furthermore plot the 

^This branch of non-uniform string solutions was in fact already discovered by Gregory and Laflamme 
in El. 

2 Other recent and related work includes PHITKIITHIITTIITHHT^I^OT) . 
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non-uniform string branch obtained numerically by Wiseman for d = 5 in the (M, n) 
phase diagram. 

A new Smarr formula for static and neutral black holes/strings on cylinders 

rg= '^~^;^ M, (1.1) 

a — 1 

is derived in Section \^ We furthermore discuss how this relation can be used to study the 
thermodynamics of branches of solutions in the (M, n) phase diagram. In particular, one 
of the consequences of the Smarr formula is an intersection rule for branches of solutions, 
stating which branch has the highest entropy. 

In Section [21 we expand on the virtues of the (M, n) phase diagram. We also compare 
to another type of phase diagram, the (A, M) diagram, that has been used in the literature 
so far. In particular, some important differences are pointed out. 

We have the discussion and conclusions in Sectional Among other things, we discuss 
the uniqueness of black hole/string solutions on a cylinder. 

A number of appendices is included providing some more details and derivations. In 
Appendix we derive the Green function for the Laplacian on a cylinder M'^^^ x S^. 
In Appendix^ we list for completeness the Gregory-Laflamme masses for uniform black 
strings in 5 < D < 10 space-time dimensions. Appendix [0 provides various details of 
the numerically obtained non- uniform solution of Wiseman [3] , which we used to plot this 
branch in our new phase diagram. 

Note added 

After this paper originally was submitted the paper |^ appeared which overlaps with 
some of the considerations of this paper. 

2 Definition of mass M and relative binding energy n on 
cylinders 

It is well known that for static and neutral mass distributions in flat space M"^ the leading 
correction to the metric at infinity is given by the mass. In this section we show that on a 
cylinder W^~^ x we instead need two independent asymptotic quantities to characterize 
the leading correction to the metric at infinity. In the following we give two such quantities, 
one of which is the mass M, and show how they can be measured asymptotically on the 
cylinder. We also determine the other independent asymptotic quantities that one can 
define. 

We note that similar asymptotic quantities as the ones given here have previously been 
considered in |2S1 123 • 

We consider a static and neutral distribution of matter which is localized on a cylinder 
M X S^. We denote here the time-coordinate as t = x^, the compact direction 
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and the directions of M . The direction z is periodic with period 2nRT- 

We also define the radial coordinate r by = {x^)'^ + • • • + (x^~^)^. With this, we assume 
the energy momentum tensor 

Too = ^>, T,, = -b, (2.1) 

with all other components being zero. Here g is the mass density and b is the binding 
energy or negative pressure along the z direction. Note that the conservation of the energy 
momentum tensor then imposes that dzb = 0. We assume in the following that g and b 
are functions of r and z only, i.e. that the matter is spherically symmetric on M'^^^. 
We define then the mass M and the relative binding energy n as 

M = J d'^x p{x) , n = ^ j d'^xb{x) . (2.2) 

n is called the relative binding energy since it is the binding energy per unit mass. We 
furthermore define Newton's gravitational potential as 

^ d-2 ^2 d"^ d-2 d , , 

V^^ = 8TTGt^——g, V2 = — + — + — . 2.3 

d — I oz^ ar^ r or 

Away from the localized distribution of mass we can then write^ 

^/ N d-2 SttGn^ I / kr\ ( kz\ , 

^) = - {d-l){d-^)2.Rr ^^-^^\-Rr) '""'[-Rr ) ' ^^''^ 

with the Fourier modes 

_{ J^^'^'^dzJ^drr'^-^pir^z) ,k = 

"'-\2jr^^dzJ-drr^-^g[^)cos(^)pir,z) , k>l ^''^^ 

where 

a — 1 \ d-3 d-3 , , , , ^ 2 2 



5(x) =r(^^j2— Id_3(x) , h{x) = ^^^^ x—Kd^{x) . (2.6) 

Here Is{x) and Ks{x) are the modified Bessel functions of the second kind. Note that 
g{0) = h{0) = 1. See Appendix^ for the Green function of the Laplacian on a cylinder 
M'*^^ X from which the above can be derived. 

We see that the Fourier modes g^ characterize the distribution of mass, and that in 
particular = M/Qd-2- Since dzb = we see that b{r,z) is completely characterized by 
M and n. Therefore, if we know M, n and g^, k > 1, we know the energy momentum 
tensor completely (when we are away from the mass distribution) . In accordance with this 
we show in the following that M, n and gk precisely correspond to the gauge invariant 
information for a given metric. 



^Notice that we assume p{r,—z) = p(r,z). This we restrict ourselves to since it appUes to all cases 
considered in Section |21 
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The Einstein equations in a {d + l)-dimensional space-time are 

R/iu -^gtj.uR = 8TrG^Ti_,„ , fi,u = 0,1, ...,d . (2.7) 
We consider now these equations for a weak gravitational field 

g/iu = rjiJ-u + hfj.u , \hfj,u\ <.l , (2.8) 

where rj^^, = diag(— 1, 1, 1) is the Minkowski metric and is the correction to the 
Minkowski metric.'^ One can consider as the first order correction to the metric 
when making an expansion of the metric in powers of Newton's constant Gn- For use in 
the following, note that our assumption that the mass distribution is static means that 
dtOfMu = and goi = 0, i = 1, d and hence Roi = 0. 
We now impose the harmonic gauge for h^u 

= iMv'^hp^) . (2.9) 
With this condition, the Ricci tensor is to first order in Gn given by 

R^.u = -\y^h^, . (2.10) 
The Einstein equations (|2.7() thus becomes 

V^/^oo = {{d - 2)g - b) , V'h^^ = {g-id- 2)b) , (2.11) 

Define now the potential B for the binding energy by 

V^B = -^b. (2.12) 
Away from the localized distribution of mass (where b{x) is also zero) we have then 



1 SvtGn nM 

(d - l)(d-3) 17d_227r i?T 



Bir) = .... ..r. .1\ -T, . (2.13) 



The unique solution of the Einstein equations (|2.11l) is therefore 

/loo = -2^ -25, h,, = -^^'^-2{d-2)B , hij = 6^, ^^<^ + 2B^ . (2.14) 

We can thus conclude that is completely determined in the harmonic gauge by M, n 
and Qk, k > 1. Obviously this means then that there are no additional gauge invariant 
parameters than M, n and g^. 



*We define h^"' = Tjt"'r]'"' h^^ so that g"" = ry""^ - /i"". 
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We now consider how to measure M, n and gi^. when we do not assume h^^i, to be in 
the harmonic gauge. To analyze this one can try to make coordinate transformations that 
take us into other gauges. However, we first need to define what coordinate transformation 
we ahow, i.e. which general boundary conditions we impose. 

Consider now the (r, z) coordinate system along with the first order correction /i^j^ to 
the metric without imposing a specific gauge. We first impose the condition that when 
Gn the zeroth order metric should be —dt'^ + dz'^ + dr'^ + r'^dQ,'^_2. Moreover, when 
r — > cxo the metric should also reduce to —dt'^+dz'^+dr'^+r'^dd'^_2- We impose furthermore 
that z is a periodic coordinate with period 2TrRx- Finally, we impose that the leading 
part of h^i, for r — > oo is independent of z. This is physically reasonable to impose since 
very far away from the mass distribution the true physical dependence on the periodic 
direction z should diminish just like it does in 1)2.4(1 where the leading term is proportional 
to Clearly this condition is satisfied in the harmonic gauge. 

Thus, if we consider an allowed coordinate transformation from {r,z) to {f,z), we see 
that it has to be of order Gn (or higher), that f/r ^ 1 and z/z ^ 1 for r ^ oo, and that 
z should be periodic with period 2ttR't- Finally, the leading part of the transformed hf^i, 
for f ^ CO should be independent of z. 

That any allowed coordinate transformation should be of order Gn, immediately has 
the consequence that /iqo is invariant under coordinate transformations, i.e. gauge invari- 
ant. 

We have in addition one more gauge invariant piece in the metric. Since h^u is indepen- 
dent of z we can easily compute that Rzz = — \'^'^hzz- Therefore, in any choice of gauge 
we have that the leading contribution to h^z for large r is of the form hzz = Cz/r'^~^ + • • • . 
If one wishes to change the constant Cz by a coordinate transformation, one needs to de- 
fine a new coordinate as z = z(l + u/r'^~^) + • • • , which gives h^z = (c^ — 2u) jr'^^'^ + • • • . 
However, clearly the coordinate z would not be a periodic coordinate with constant pe- 
riod IttRt anymore so this is not an allowed coordinate transformation. Therefore, Cz is 
invariant under the allowed coordinate transformations. 

In conclusion, we have shown that 



independently of the choice of gauge. Using this we can measure M, n and Qk, k > 1, 
from a given metric. 

In particular, if one wishes to measure the mass M and the relative binding energy n 
from a given metric which has the leading behavior 




(2.15) 



to first order in Gn, and that the leading correction to Qzz for r ^ oo is 





goo = -1 + 



9zz = 



1 + 



(2.17) 
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for r — > oo, then M and n are given by 

These formulas enables us to read off M and n for any static and neutral solution. 
Note that from M and n we can define the tension in the circle direction as 

/i = — , (2.19) 

where L = 2'kRt is the circumference of the circle. The tension ^ has previously been 
considered in pHl OHl OH] . 

Bounds on the relative binding energy 

We have the following general bound on the relative binding energy 

0<n<d-2, (2.20) 

which holds for any static and neutral solution. 

The bound n > was found in j25ll26j . The proof is similar to the simple proof of the 
Positive Energy Theorem of |27j which was generalized to include horizons in |2H| • One of 
the basic assumptions that goes into the proof is the Dominant Energy Condition. 

The bound n < d — 2 follows instead from the Strong Energy Condition. The Strong 
Energy Condition gives that Too + ^^7"'^ > which in the above analysis of Newtonian 
matter gives q > from which n < d — 2 follows. Basically the Strong Energy Condition 
imposes in this case that gravity is not repulsive asymptotically. 



3 M and n as variables for phase diagram 

In this section we consider static and neutral black holes and strings on a cylinder. We 
have shown above that for any neutral and static object on a cylinder M^^^ x we can 
measure the mass M and the relative binding energy n. It is therefore natural to use 
these two independent quantities as the variables in a two-dimensional phase diagram for 
static and neutral black holes and strings on a cylinder. As we review below we need a 
two-dimensional phase diagram to illustrate the known solutions. 

We have three known branches of solutions corresponding to neutral and static black 
holes and strings on a cylinder. First we have the black hole branch. This branch starts at 
(M, n) = (0, 0), since it is physically clear that a small mass black hole will have negligible 
relative binding energy (see |S] for quantitative discussion of this) . The black hole branch 
should terminate at a certain critical mass since the horizon radius at some point becomes 
too large to fit on the cylinder. This branch is the subject of |H], where an ansatz for 
solutions in this branch was given. We have made a qualitative illustration of this branch 
in the phase diagram in Figure E 
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The second branch is given by the uniform neutral black string solutions, which exist 
for any mass M. These solutions are obtained as a direct product of a Schwarzschild black 
hole solution with a circle, i.e. the metric is 

d-3 



-fdt" + r^dr^ + + r2dl72_^ / = 1 - ( ^ 



(3.1) 



From the definition H2.17() we infer that = for any M, so that it follows from ()2.18|) 
that n = l/{d — 2) for all M. We also illustrated this in the phase diagram in Figuren(we 
are considering the special case d = 5 in the diagram). As was discovered by Gregory and 
Laflamme in the solutions in this branch are classically stable when M > Mql and 
classically unstable when M < Mql, where Mql is found numerically in pj (in Appendix 
^we have listed the Gregory-Laflamme masses for d = 4, ..,9). 

The third branch is the one discovered by 15 and further explored in |TJ . This 
branch of solutions begins with the uniform black string solution in (M, n) = (Mgl, ^/{d— 
2)) and continues with increasing M and n (at least for the known cases d = 4, 5). We have 
incorporated this in the phase diagram in Figure ^ for the case d = 5 by using the data 
points found by Wiseman in [3]. We refer to Appendix O for details on this numerically 
obtained solution. 

n 



0.3 



0.2 



0.1 



Uniform string branch 



Non-uniform string branch 



Biocl< hole branch 







2M, 



M 



Figure 1: (M, n) phase diagram for d = 5 containing the uniform string branch and the 
non-uniform string branch of Wiseman. The black hole branch is sketched. 



We thus see that all the known solutions corresponding to black holes or black strings 
can be mapped in our (M, n) phase diagram, as done in the phase diagram of Figure ^ 



4 Thermodynamics in the (M, n) phase diagram 

In this section we first derive a new Smarr formula for static and neutral black holes/strings 
on cylinders. We then show how one can use this to study the thermodynamics of branches 
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of solutions in the (M, n) phase diagram. 

To derive the new Smarr formula, we begin by considering the Komar integral 

Is = -T7^ [ dS^.D^'k^ , (4.1) 



where 5" is a d — 1 dimensional hypersurface in the d + 1 dimensional cylinder space-time 
and fc is a killing vector for the metric. We can also write this integral as 

= - Ifsk L ^(d^W A - A </x-. . (4.2) 

Consider now a static neutral solution on the cylinder with an event horizon. Consider 
furthermore a certain time t = to. Define Sh to be the null-surface of the event horizon at 
t = tQ. We also choose a d — 1 dimensional surface at r = oo for t = to, which we call Sq^, 
so that essentially Soo = 5"^^^ x S^. By Gauss theorem we have 

where V is the d dimensional volume between Sh and Soo so that dV = ShU Soo- Since 
for a killing vector we have D^D^k'^ = R'^''uk^ and since = everywhere in V , as 
we are away from the black hole/string, it follows that D^D^k'^ = 0. Thus, we find that 
= Is^ ■ 

Since we have a static solution we can choose k to be the time-translation killing 
vector field, i.e. k = d/dt. On the event horizon Sh we can choose a normal n such that 
kfj.n^^ = —1. Then one can show that dS^y = {k^riy — ki,ny)dA where dA is the area 
element on the event horizon (see |29ii). Using this, we have j29| 

where we used that k^Dj^k^ = nk^. Here k is the surface gravity, k^n^ = — 1 and we used 
the thermodynamic identifications k = 2'kT and A = AG^S. On the other hand, we have 
asymptotically 

where we first used that the non-zero components of D^k^ are D^k^ = —D^k^ = ^(9r5oo 
to leading order, and (|2.18j) in the last step. 

In conclusion, we get from ()4.4|) and (|4.5() that 

TS = 'tLl^M . (4.6) 
a — 1 

This is a new Smarr formula which holds for any static and neutral black hole/string on 
a cylinder. Indeed the Smarr formula (|4.6() can be checked for all three known branches 
in the phase diagram. For uniform black strings {n = 1/ {d — 2)) as well as for small black 
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holes on the cyhnder (n = 0) it reduces to known Smarr relations. In Appendix [0 we 
show that the formula also holds to very high precision for the non-uniform black string 
branch using the numerical data of Wiseman. 

Define now a branch to be a series of solutions corresponding to a curve in the (M, n) 
diagram for which the first law of thermodynamics 5M = T5S is obeyed when moving 
between the solutions.^ Then if we consider a branch for which we know the precise 
curve in the (M, n) phase diagram, we can find the thermodynamics on the branch. For 
simplicity we consider here the case where the curve is given by a function n = n{M). 
Using now 5M = T5S and the Smarr formula (|4.()j) we get 

6\ogS ^ 1 d-l 
SM Md-2-n{M) ' ^ ' ' 

Thus, we see we only need to know the entropy in one point of the branch to find the 
complete function S{M) for the branch.^ A similar statement applies to the temperature. 
In conclusion, knowing the exact curve of a branch in the (M, n) phase diagram enables 
one to obtain the entire thermodynamics of that branch. 

Another useful feature of depicting a branch in the (M, n) phase diagram is that when 
two branches intersect we have the following Intersection Rule stating which branch has 
the highest entropy.^ The proof of the Intersection Rule is at the end of this section. 

< Intersection Rule 

Consider two branches 71.2 in the (M, n) phase diagram, parameterized by the func- 
tions ni(M) and n2(M). Let the range of masses be from M' to M" (M' < M") 
and assume that ni(M) > n2(M) for M' < M < M" . Then 

71 intersect 72 at M = M' Si{M) > 52 (M) for M > M' , (4.8) 

and 

71 intersect 72 at M = M" ^ Si{M) < S2{M) for M < M" . (4.9) 

Thus, we see that for two intersecting branches we have the property that for masses 
below an intersection point the branch with the lower relative binding energy has the 
highest entropy, whereas for masses above an intersection point it is the branch with the 
higher relative binding energy that has the highest entropy. We have illustrated this in 
Figure 

^We show in I3UI that all series of solutions which are connected by infinitesimal steps (i.e. being 
differentiable) have the property that the first law SM = TSS is always obeyed. In other words, for a 
smooth curve in the (M, n) phase diagram the first law is obeyed on the curve. 

"Notice however the important subtlety that if the point in which one knows the entropy has zero 
entropy then more information is needed to get the entropy function S{M) since it is evident from 14.711 
that one can only get log S{AI) up to an additive constant in that case. 

^Notice however that if the branches intersect in a point for which the entropy is zero the Intersection 
Rule is not necessarily true. If we denote M* the mass corresponding to the intersection then we need that 
Si{M)/S2{M) 1 for M M* in order for the Intersection Rule to hold. 
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-► 

M 

Figure 2: Illustration of the Intersection Rule. 



As a consequence of the Intersection Rule, we find that it is impossible for two branches 
which are described by well-defined functions n = n{M) to intersect twice. If we write 
the two functions corresponding to the two branches as ni(M) and n2(M) and the two 
intersections as M' and M", then the Intersection Rule gives a clear contradiction since 
we get from (gSl that 5i(M") > S2{M") and from ^ that 5i(M') < S2{M') both of 
which contradict the fact that the two branches intersect in M' and M" . 

We remind the reader of the physical reason why we are comparing the entropy of two 
branches at a given value of M, i.e. in the micro canonical ensemble. This is because, in 
order to address the stability of a particular branch, one would like to examine whether 
it is entropically favorable for the solution to move to another branch, characterized by a 
different value of the relative binding energy, by redistributing its mass. 

The proof of above stated Intersection Rule goes as follows. Let 5j(M) = S{M,ni{M)), 
i = 1, 2 be the entropy on the two branches 71^2- If the two branches meet at M = we 
can write the relation 



52 (M) Jm, 6M \S2{M) 

Obviously is either M' or M" . Consider now the differential 



5M\S2{M)J T1T2SI ' '{d-l)TiT2Sl 

where we used the first law of thermodynamics 5Mi = Ti5Si for each of the branches 
in the first step and the Smarr formula (|4.6|) in the second step. Since we assume that 
ni(M) > n2(M) for M' < M < M" we have 

^ ( > for M' < M < M" . (4.12) 

5M \S2{M) J ^ ' 

Using (|4.12|) along with (|4.1fl|l . we see that putting = M' gives (|4.8|) and putting 
= M" gives (|4.9() . completing the proof. 
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5 Virtues of the (M, n) diagram and comparison to the (A, M) 
diagram 



So far in the literature, the phase diagram used to display the various phases of strings 
and black holes on cylinders has been what we call the (A, M) diagram. It is defined using 
the so-called non-uniformity parameter of [1] 

A = ^f|^-l), (5.1) 

where iimax and i?min are defined respectively as the maximal and the minimal radius of the 
S'^~'^ sphere on the horizon. Note that A is invariant under coordinate transformations.^ 
In the following we discuss some important differences between the (A, M) and the (M, n) 
phase diagrams. 

An immediate physical point to notice is that M and n completely characterize the 
leading correction to the space-time at infinity. Thus, if an observer at infinity should 
draw diagrams of his/her findings, it would be natural to draw these in an (M, n) phase 
diagram. In contrast, A can not be measured at infinity, since to know A one has to know 
the space-time all the way to the horizon. 

Note also that the parameter n always stays finite contrary to A which is infinite for 
the black hole branch since i?min = on this branch. Moreover, we have already noted in 
Section [21 that < n < d — 2. However, we do not expect n > l/{d — 2) since it seems 
reasonable to assume that the uniform string branch, which has n = l/{d — 2), has the 
maximal possible relative binding energy. We therefore have the bound < n < ^j^. 
That n is bounded makes it useful as a measure for "how far away" a given string solution 
is from a given black hole solution. 

However, we should note that the (A, M) diagram has information that the (M, n) 

diagram does not: The value of A explicitly states whether a solution is a uniform string, 

a non-uniform string or a black hole. On the other hand, as we described in Section |3 

the (Af , n) phase diagram has several useful features. For example we have explained how 

one can read off the thermodynamics of a branch of solutions using the (M, n) diagram. 

Moreover, it was explained that there are certain rules in the (M, n) diagram restricting 

which branches are possible. 

*This can easily be seen by considering the most general static metric with spherical symmetry in the 
R''"^ part. Write here this metric as 

ds^ — ~Udt^ + Qijdx^dx-' + VdVl^ci-2 , 

with i,j = 1, 2 and with U, V and gtj being functions only of and x^. Then the size of the S''~^ sphere 
at a point on the horizon is given by the value of V at that point. It is then clear that the maximum and 
minimum of V is invariant under any coordinate transformation {x^,x^) (x^,i^) since V transform as 
a scalar. 
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6 Discussion and conclusions 



We have discussed in this paper a novel type of phase diagram for black holes and black 
strings on cylinders. The phase diagram involves, beyond the mass M, a new asymptotic 
quantity n, which is the relative binding energy in the circle direction. So far, three 
branches of solutions are known in the phase diagram: the uniform black string, the non- 
uniform string solutions numerically obtained by Wiseman [SI (see also jl]) and the small 
black hole branch. We have plotted these known branches in the new phase diagram. 

We have derived a new Smarr formula, valid for black holes and strings on the cylinder, 
which involves the quantity n. This Smarr formula provides a powerful instrument as it 
enables one to obtain the entire thermodynamics, once the solution n{M) in the phase 
diagram is known. Moreover, when combined with the first law of thermodynamics, the 
Smarr formula also implies an Intersection Rule, which easily determines the branch of 
highest entropy for intersecting branches in the phase diagram. Although we have fo- 
cussed in this paper on neutral black holes /strings, this part of the development is easily 
generalized to the non-extremal and near-extremal case as well as the inclusion of angular 
momenta. For example, in the non-extremal case, the Smarr formula is simply obtained 
by M ^ M - /iQ in 

We also remark that the formalism introduced in this paper seems well-suited for 
the important study of the phase structure of black p-branes in spaces with more than 
one compact direction. The simplest generalization one can think of would be to study 
black holes, strings and membranes on the torus W^^"^ x T^. Following our development, 
one would need in that case, beyond the mass M, three additional asymptotic quantities 
corresponding to the stress in the torus direction. 

We have pointed out some of the virtues of the (M, n) phase diagram and compared 
to the (A, M) diagram, so far used in the literature. In a forthcoming work [2111 ^'^^ 
use the formalism introduced in this paper to obtain various new insights into the phase 
structure for black holes and strings on cylinders. 

Finally, we comment on the issue of uniqueness of black hole/string solutions on a 
cylinder. We obtained in Section [2| the gauge invariant quantities that determine the 
leading order metric in the presence of a static and neutral (spherical symmetric) mass 
distribution on the cylinder. Beyond the mass M and the relative binding energy n, these 
were shown to also involve the higher Fourier modes Qk>i of the mass distribution. For 
black holes/strings on cylinders this implies the possibility of having several solutions for 
a given point in the (M, n) phase diagram, distinguished by the Qk- 

For neutral black holes in four space time dimensions it is well-known that the solu- 
tion is uniquely characterized given the mass.^ It would naturally be very interesting to 
examine whether similar uniqueness theorems hold for black holes/strings on cylinders.^'' 

®For asymptotically flat solutions in higher dimensions no simple uniqueness theorems are presently 
known. In particular, in Ref. (31| a new black ring solution was found in five dimensions. 

^'^We recall here that our analysis pertains to cylinders R'*~^ x 5^, d > 4, so that the number of space 
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For example, it may be that the Fourier modes are determined in terms of M and n. 
It is thus conceivable that the following uniqueness hypothesis is valid: 

■ Consider solutions of the Einstein equations on a cylinder M'^-i x S^. For a given 
mass M and relative binding energy n there exists at most one neutral and static 
solution of the Einstein equations with an event horizon. 

A less strong version of the above hypothesis would be that also the choice of horizon 
topology (S"^~^ X or S^^^) needs to be specified in order to obtain a unique solution. 
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A Green function on a cylinder 

We consider the cylinder M^~^ x with metric 

ds'^ = dz^ + dr^ + r^dnl_2 . (A.l) 

where z is periodic with period 27r Rt- Here the S'^~'^ sphere is parameterized by the angles 
9i, ...,9d-2 and it has the metric 

d-2 

d^l-2= Yl ^abdeadOb . (A.2) 

a,b=l 



Define 7 = y^det(7ab). Define furthermore the notation Q = {61, ...9(i_2)- 
We define then the Green function on M'^-i X S as 

Vl^^eG{z -z',r- r', &-&') = -^^6{z - z')6{r - r')6{e - &') . (A.3) 

In the following we derive a useful expansion of the Green function on the cylinder which 
is used in Section [21 

We first introduce the spherical harmonics Yi{Q) which are eigenfunctions 

d-2 

-Y^dea \ll''''de,Wu-^^d-2) \ = -hiYi{e,, ...,9d-2) , (A.4) 



time dimensions is five or higher. 

^^Note that in [ISJ a difi'erent proposal for a uniqueness theorem was given: For a given mass AI there 
exists only one stable solution of a certain topology. That proposal is of a rather different nature, since 
if correct it does not address the point of having different branches that exists classically. Whether this 
proposal is correct or not is thus independent of the above discussion of the uniqueness properties relating 
to Ai and n. 
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that satisfy the completeness relation 



(A.5) 



We can now expand any function on S'^ ^ as f{Q) = Xl/^'^K©)- We can furthermore 
expand the delta function for z as 



2-kRt 



cos 



fe=l 



k{z - z') 



The Green function can then be expanded as 

G{z-z\r-r',Q-Q') = ^ (6')^/ (6) m/'^^ (r, /)e^^(---')/i?T 



Using this the Green function equation ()A.3|) then becomes equivalent to 



(A.6) 



(A.7) 



d-2^ \ . nAk)(^ - -^^6{r - r') . (A. 



In order to solve this equation, we notice that the modified Bessel functions of the second 
kind Ky{x) and Ii,{x) are solutions to the equation 

,2^ 



Notice then that if we define g as 

g{x) = x~^f{kx) , 
then (7 is a solution to the equation 



9 H 9 

X 



So, using 



1^1 



-{d-sy + hi 



9 = 0. 



(A. 



(A.IO) 



(A.ll) 



(A.12) 



we see that g{r/RT) is a solution of l|A.8|) when r ^ r' . Define now the two independent 
solutions of the differential equation l|A.ll|) by 



ipi {v^x) = k ^x 2 luikx) , -02 (^;^) = k'^^ ^ Ki,{kx) , 



(A.13) 



satisfying 



iIji^^\v,x) ~ 



1 



x""-'-^ , V'fV'^) - 2''"^r(i/)x"''""2" , x<.l (A.14) 



2''T{v + l] 



(A.15) 
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Then we can write the solution of HA.8|) as 

^d^2 j(k) 



Ml^\r, r') 



nd~-3 

Itrp 



with 




(A.16) 



(A.17) 



B Critical masses for the Gregory- Laflamme instability 

For completeness we list in this appendix the Gregory-Laflamme masses for uniform black 
strings on a cylinder W^~^ x with d = 4, 9. The data has been taken from (the 
published version) by reading off /i from their Fig. 1 . The values of that we read off are 
0.44, 0.62, 0.80, 0.93, 1.04, 1.15 for d = 4, 9. One can then see in P H] that the critical 
horizon radius is given by vq/Rt = 2;U, where ro is the horizon radius of the uniform black 
string solution H3.1|) and Rt is the radius of the cylinder. Using that 



M 



nd-227rRT 



(d-2)rl 



(B.l) 



levrGAT 

we can compute the critical masses. In Tabled we have listed the horizon radii and the 
masses for d = 4, ...,9. We have used the notation that the circumference of the cylinder 
is L = 27ri?r. 



d 


4 


5 


6 


7 


8 


9 




0.14 


0.20 


0.25 


0.30 


0.33 


0.37 


167:GnMgl/L'^~^ 


3.52 


2.31 


1.74 


1.19 


0.79 


0.55 



Table 1: Critical horizon radii and masses for the Gregory-Laflamme instability. 



C Details on Wiseman's non-uniform black string data 

In this appendix we provide the reader with details of how the plot of Wiseman's non- 
uniform black string branch in Figure ^ was made. Wiseman's non-uniform black string 
solutions (with d = 5) were obtained in ^ by numerically solving Einstein's equations for 
the conformal ansatz 

dsl = ^e^^dt^ + e2B/^^2 ^ ^^2. ^ g2C. ^ r^)dn^ . (C.l) 

Here, A,B,C are functions of (r, z), with asymptotic forms 

A^^, Bc.^, + (C.2) 
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Figure 3: Plot of 02 and 62 versus A using the data of 0. 



The data found in 5 consist of 25 points, each corresponding to a non-uniform black 
string solution. Each of the 25 points is then specified by six numbers: A, 02, 62, T-^, 
and M^. A is defined in (|5.1() . T^, and M„ are the temperature, entropy and mass 
that was computed in 5 , and 02 and 62 are defined in (1^121) . These dat a was generously 
made available to us by Toby Wiseman. 




Figure 4: Plot of M and n versus A. AI and n are computed from the 02 and 62 data. The crosses 
mark the actual data-points. 

The only data necessary to make the phase diagram in Figure ^ are 02 and 62 for the 
25 points. We have plotted 02 and 62 versus A in Figure |31 We refer the reader to ^ for 
plots of T„, S'w and M„ versus lambda. By comparing 1)2. 17(1 and (|C.1I) - ()C.2|) we see that 
a2 and 62 are related to ct and as 



ct = 1 - 2a2 , Cz = 2b2 , 



(C.3) 



where we note that m = 1 in ()C.1|) . We can then use this in (|2.18() (with d = 5) to 
determine the mass M and relative binding energy n for each data point according to 

M ct-3cz . . 

Mgl 3 3q - 

In Figure 13 we have plotted M and n versus A. Moreover, in Figure [S] we have plotted M 
versus n. 

We now check that the first law of thermodynamics 6M = T5S is obeyed for the above 
computed M. In left graph of Figure El we have depicted 6M/5X and T6S/6X versus A. 
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We note that they agree to a high precision. The third curve in Figure El is 5Mw/5\ 
where Myj is the mass that Wiseman computed in the original version of [Sj. We see that 
the systematic error that Wiseman noted in the original version of |51 in fact is due to 
an erroneous computation of the mass, and not an inaccuracy of the solutions that were 
found (see also comment in new version of [5J. 




1 2 3 AO 1 2 3 4X 



Figure 6: On the left plot we have checked the first law of thermodynamics for Wiseman's branch 
of solutions. Here the solid curve shows the crosses show ^ and the diamonds show ^|^, 

all in units of Mql- On the right plot we have checked Smarr's formula (|4.6() . Here the solid curve 
show TS while the crosses show ^^M, both in units of Tql'S'gl = f^'^GL- 

In the right graph of Figure IHl we have instead plotted TS (computed from and 
S^) and versus A. This is a check of Smarr's formula 1)4. 6|) . We see that Smarr's 

formula is correct to a very high precision. 
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